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Abstract. We analyze the free energy and construct the Gibbs-KMS states for a class 
of quantum lattice systems, at low temperature and when the interactions are almost 
diagonal, in a suitable basis. The models wc study may have continuous symmetries, 
our results however apply to intermediate temperatures where discrete symmetries are 
broken but continuous symmetries are not. Our results are based on quantum Pirogov- 
Sinai theory and a combination of high and low temperature expansions. 
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1. Introduction 

In this paper we study the low temperature phase diagram for a class of quantum 
lattice systems. Starting with |PS[ |Sin|| , Pirogov-Sinai theory has evolved [ KP , ^ah , BKL , 



BS , Bl , BK | into a very powerful tool to study the pure phases, their coexistence and 
the first-order phase transitions in classical spin systems at low temperature. In recent 



years large part of the Pirogov-Sinai theory has been extended to quantum systems [Pir, 
BKU, DFF, DFFR, [KU|], quantum spin systems as well as fermionic and bosonic lattice 



gases, and applied to a variety of models |FR, DFF2 , GKU ] to describe insulating phases 
associated with discrete symmetry breaking. Here we formulate the Pirogov-Sinai theory 
in terms of tangent functionals to the free energy. This allows to discuss the completeness 
of the phase diagram avoiding the difficulties associated with boundary conditions. We 



reformulate results of [BKU, DFF, DFFR, KU] in this framework, and extend the theory to 
a class of models where discrete symmetries are broken at intermediate temperatures. This 
applies in particular to some systems with continuous symmetries. For this, we consider 



the restricted ensembles introduced in [ BKL ] that are very useful to analyze phases which 
are associated to a family of configurations rather than to a single configuration. 
The models that we consider have Hamiltonians, for finite volumes A, of the form 

H A = V A + T A 

where V\ is a classical Hamiltonian (i.e. diagonal in a suitable basis) and T A is a (usu- 
ally small) quantum perturbation. In typical situations the suitable basis is the basis 
of occupation numbers of position operators. Electronic systems provide a large class of 
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interesting models. The classical interaction V\ describes the many-body short range and 
classical interaction between the spin-^ fermions as well as external fields and chemical 
potentials: 

Va — ^ ~] Jx,a^x,a ^ ^ J xy^aa'^x^a^y^a 1 ~\~ ' ' ' • 

x£A x,y£A 

A typical quantum perturbation T\ is the kinetic energy 

— ^ ^ txy,a(Cx<jCy<j ~\~ h.C.) , 
<X,y>CA 
<r€{U} 

where c X(T and c xa are the creation and annihilation operators and <x,y> denotes pairs 
of nearest neighbors. 

Often, in such systems, the behavior at low temperatures arises from a subtle interplay 
between the (classical) potential energy and the kinetic energy. In this paper two such 
mechanisms are considered and combined, each of which we now illustrate with an example. 

Example 1: Hubbard Model. In this case the (classical) interaction is only on-site: 

Va = y^ j Un x ^n xl - /i(n x f +n xl ). 
x£A 

For suitable values of U and //, the ground states of V\ have an infinite degeneracy (in the 
thermodynamic limit): each site is occupied by a single particle of arbitrary spin. However 
the kinetic energy lifts this degeneracy and induces an effective antiferromagnetic interac- 



tion between nearest neighbors. The perturbative methods of [DFFR, DFF2] shows that, 
in this parameter range, this system is equivalent, in the sense of statistical mechanics, to 
the Heisenberg antiferromagnet, up to controlled error terms. If the hopping coefficients 
are asymmetric (e.g. t xy t <C t xy \ ) then quantum Pirogov-Sinai implies the coexistence of 
two antiferromagnetic phases at low enoug h temperatures ||DFFR| , |KU|, pFF2| . Ri gorous 
results for the Hubbard model are reviewed in [Lieb]. 



Example 2: Extended Hubbard Model. This variant of the Hubbard model includes 
a nearest neighbor interaction: 



V A = ^2 I Un x1 n xl ~ At(?Vr + n>xl) \ +W (n x1 + n xl )(n y1 + n y[ ) . 

xeA ' <x,y>cA 

If the the interaction between nearest neighbors is repulsive then for suitable values of 
U, W and /U the ground states of V\ are chessboard configurations where empty sites 
alternate with sites occupied with one particle of arbitrary spin. The degeneracy of the 
ground states is infinite in the thermodynamic limit but we have a spatial ordering of the 
particle. Using restricted ensemble we associate a pure phase to this spatial ordering by 
neglecting the spin degrees of freedom. The methods of this paper imply the existence of 
only two pure phases in the intermediate temperature range 

(3t < 1 and /3W > 1 . 

The temperature is so low that the spatial ordering of the particles survives but so high 
that the spin are in a disordered phase. The continuous symmetry (if t xy ^ = t xy ^) is not 
broken in this parameter regime. 
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Figure 1. Spiral order in 



These two models illustrate some of the mechanisms arising from the competition be- 
tween classical and quantum effects, where the system remains insulating and no con- 
tinuous symmetry is broken. Our main result, Theorem 4.4, provides tools to describe 



the phase diagram of such models, in particular the coexistence of several phases and the 
associated first-order phase transitions. 

The main technical ingredient in this paper is a combined low-temperature and high- 
temperature expansion for suitable contour models obtained using the perturbation theory 
developed in ||DFFR| . 

This paper is organized as follows. In Section || we describe the general formalism of 
quantum lattice systems and the perturbation theory of [ DFFR |. Section || is devoted 
to the Pirogov-Sinai theory. In Section |] we state the results of Pirogov-Sinai theory for 
quantum systems. The extended Hubbard model is discussed in Section || as an illustration. 
In Section ^ we prove our main result by studying a contour model and deriving the 
required bounds on the contours. 

2. General framework of quantum lattice models 
2.1. Basic set-up. We consider a quantum mechanical system on a i/-dimensional lattice 



Z", as considered, e.g., in [ Rue , Isr , BR , Sim | . We will need a slight modification of the 
usual formalism in order to treat fermionic lattice gases [DFFR] and to accommodate the 
fact that fermionic creation and annihilation operators do not commute but anticommute. 
A quantum lattice system is defined by the following data: 

(i) Hilbert space. For convenience we choose a total ordering (denoted by the symbol 
X) of the sites in We choose the spiral order, depicted in Figure |l] for v = 2, and an 
analogous ordering for v > 3. This ordering has the property that, for any finite set A, 
the set A := {z £ Z u , z ^ A} of lattice sites which are smaller than A, or belong to A, is 
finite. To each lattice site a £ Z y is associated a finite-dimensional Hilbert space Tt a and, 
for any finite subset A = {a± ~< • • • -< a n } C 7L V , the corresponding Hilbert space Ha is 
given by the ordered tensor product 



Ha = H 



We further require that there be a Hilbert space isomorphism <j) a : H a 



(2.1) 
H, for all 



(ii) Field and observable algebras. For any finite subset A C W an operator algebra 
Ta, the field algebra, is given. The algebra Ta is isomorphic to the algebra B(Ha) of 
bounded operators on Ha, but in general Ta ^ B(TCa), rather Ta C B(H-^). The algebra 
J- A is a *-algebra equipped with a C*-norm obtained from the operator norm on B(Ha)- 
If A C B and a -< b, for all a G A and all b E B \ A, then there is a natural embedding of 
Ta into Tb- An operator K G Ta corresponds to the operator K g5 1h b \ a m J~B- In the 
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following we denote by K both operators. For the infinite system the field algebra is the 
C* algebra given by 



F = U ^ , (2.2) 

A/Jf 

(the limit being taken through a sequence of increasing subsets of Z", where increasing 
refers to the (spiral) ordering defined above). 

The algebras J- a contain the observable algebras Oa which have the same embedding prop- 
erties as the field algebras and, moreover, satisfy the following commutativity condition: 
If A n B = 0, then for any K G T A , L G Ob we have 

[K,L}=0. (2.3) 

For the infinite system the observable algebra O is given by 



O = |J O a - (2.4) 

A/"E V 

The group of space translations 7L V acts as a *-automorphism group {T a } a &'L v on the 
algebras T and O, with 

F x (Fx), O x (Ox), (2.5) 

for any X C Z" and a G Z". 

(iii) Interactions, dynamics and free energy An interaction _£f = {Ha} is given: 
This is a map from the finite sets A<z7L v io self-adjoint operators Ha in the observable 
algebra Oa- We assume the interaction to be translation invariant or periodic, i.e., there 
is a lattice T C Z", with dimT = z/, such that r a HA = Ha+o,, for all a G T and all 
A C Z". We will consider finite range or exponentially decaying interactions. The norm 
of an interaction is defined as 

||if|| r = sup V \\H A \\e rlAl , (2.6) 

a ^Ata 

for some r > 0. Here \A\ denotes the cardinality of the smallest connected subset of 7L V 
which contains A. We shall denote by B r = {H : ||JjT||f < oo} the corresponding Banach 
space of interactions. 

For a finite box A, we denote H\ the finite-volume Hamiltonian given by H\ = 
Y2acA Ha- Here, we consider only periodic boundary conditions, i.e. A is the j/-dimensional 
torus (Z/XZ)^, L being the size of A. In the sequel we will consider infinite volume limits; 
the notation lim\y^ will stand for lim/^oo. 

If H G B r , the interaction H determines a one-parameter group of *-automorphisms, 
{at}j g iR on T . These automorphisms are constructed as the limit (in the strong topology) 
of the automorphisms given by for K G Fa, A C A by 

a^(K) = j tH *Ke~ itHA . (2.7) 



The proof is standard (see e.g. yBR| j). Note that one makes crucial use of the commuta- 
tivity condition Q2.3] ). 

For an interaction H and at inverse temperature (3 the partition function is defined as 

ZjJ = Tr e _/3HA ; (2.8) 
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the free energy f(H) is then 



Existence of the limit is a well-known result, see |Isr| , |Shn]. Notice that f(H) is a concave 
function of the interaction H. 

(iv) KMS states and tangent functionals. A state w on O is a positive normalized 
linear functional on O. A state w is periodic if w o r a = w, for all a in a lattice V and 
invariant if T = Z". A KMS state at inverse temperature is a state wp which satisfy the 
KMS condition 

wp(Ka t {L)) = w p {a t .^{K)L) . (2.10) 

For finite systems with periodic boundary conditions it is easy to check that the Gibbs 
state given by 

w pA (-) = (Tre-^A)-i Tr(e -^A .) (2 .ii) 

satisfies the KMS condition. The set of KMS states is convex, and w is called extremal if 
it cannot be written as linear combination of KMS states. The state w is clustering if 

lim w(K Ta (L)) = w(K)w(T a L) , (2.12) 

a— too 

for all K, L E O. Note that a state w is extremal if it is clustering. The state w 
is exponentially clustering if, for any local observables K £ Oa, L G Ob we have the 
property 

w{Kr a {L)) - w{K)w{r a L) ^ C k ,l^ HK (2.13) 

with £ > 0; here Ck,l depends on K and L only. 

If we consider the free energy as a function of the interaction, KMS states at inverse 
temperature f3 are in one-to-one correspondence with tangent functionals to the free energy. 
The free energy / is a concave function of the interaction H and a linear functional a on 
B r is said to be tangent to / at H if for all interaction K £ B r we have 

f(H + K) < f(H)+a(K). (2.14) 

To an invariant state w we associate a tangent functional a defined by 

a(K) = w(A K ) (2.15) 

where Ape = ^xso l^l -1 -^^ (and similarly for periodic states). The results of Israel and 
Araki | jlsr| , Araf show that if a is a tangent functional at H, then the invariant state w 



defined in ( fU5|) is a KMS state at temperature and, conversely, for any KMS state at 



temperature (3 there is a unique tangent functional a. The identification of KMS states 
with tangent functionals will be very useful to describe the phase diagrams arising from 
Pirogov-Sinai theory. 

Example. As an illustration of the general formalism we consider spin 1/2 fermions, as 
in the examples treated in this paper. The Hilbert space TL a is isomorphic to C 4 . We 
denote c\ a and c au the creation and annihilation operators of a particle at site a with spin 
a £ {Til}- O ne can construct an explicit representation of the creation and annihilation 
operators as operators in B(Ha), see e.g. Section 4.2 in [DFFR], but c aa ,cta ^ B(H a ) 



The algebras Ta C B(7i-j) are chosen to be the algebras generated by c acr , eta, a £ A, 
o G {1, 1}. The observable algebra Oa are chosen as the algebras generated by pairs of 
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creation or annihilation operators. It is easy to check that the elements J- a and Oa satisfy 
the commutativity condition ( |2.3| ). 

Classical interactions A particular class of interactions consists of the classical interac- 
tions. Let {ej}j£i be an orthonormal basis of H. Then, for A C Z^, 

£a = {®aeAe%} , with e a ja = ^ ej , (2-16) 

is an orthonormal basis of Ha- We denote by C(£a) the abelian subalgebra of Oa consisting 
of all operators which are diagonal in the basis £a- An interaction V is called classical, if 
there exists a basis {ej}j^j of H such that 

V A G C{£ a ), for all A C 1? . (2.17) 

The set Qa of configurations in A is defined as the set of all assignments {ja}{ a &A} °f an 
element j a G / to each a. A configuration uja is an element in VI a- There is a one-to-one 
correspondence between basis vectors ® a eA e j a °^ ^-A and configurations on A: 

(g)e| a ^w A = {ja}aeA. (2.18) 

a£A 

In the sequel we shall use the notation e^ A to denote the basis vector defined by the con- 
figuration uja via the correspondence ( [2.18 ). Since a classical interaction V only depends 



on the numbers 

$ A (u A ) = (euJl^lO (2.19) 

we may view §a as a (real-valued) function on the set of configurations. Similarly the 
algebra C(£a) may be viewed as the *-algebra of complex- valued functions on the set of 
configurations £Ia- 

2.2. Perturbation theory for interactions. The interactions we will study have the 
form H = V + AT where V is a classical interaction, T is a perturbation and A a small 
parameter. A typical situation is the following: the classical part of the interaction has 
infinitely many ground states, i.e the number of ground states of the finite-volume Hamil- 
tonian H\ diverges as |A| — > oo, but the perturbation V lifts this degeneracy (completely 
or partially). This is usually easy to check this using standard perturbation theory for the 
finite- volume Hamiltonian Va + ATa. Standard perturbation theory however does not work 
in the thermodynamic limit, the norm of the error growing with [A| and other methods 



are required. Such methods have been developed in [DFFR] and applied in [FR, DFF2] 
(see also |KU[1 for an alternative approach). 

The idea is to construct an interaction H which is equivalent to H and which can be 
cast in the form 

H = V{\) + T(A) , (2.20) 

where now the degeneracy of the ground states of V is lifted and T(A) is suitably small 
with respect to V(X). 

Recall that two interactions H and H are equivalent if there exists a *-automorphism 
of the algebra O of local observables such that 

H A = l(H A ) , (2.21) 

for all A. In particular, if H G B r , there exists f such that H G Bf- A convenient way of 
constructing equivalent interactions is with a family of unitary transformations U\. Let 
Sa, A C Z u , be a family of antis elf adjoint operators, periodic or translation invariant, with 
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S A £ Oa and \\S\\ r < oo for some r > 0. We set Sa = XmcA &a and then U\ = exp(SA) 



is unitary. It is shown in [DFFR] that if ||S|| r is small enough then the unitary equivalent 



Hamiltonians Ha = UaHU^ 1 define an interaction H G Bf for some f > and H is 
equivalent to H. 

We consider now an interaction of the form H = V + AT which satisfy the following 
conditions 

(PI) The interaction V is classical and of finite range. Moreover, we assume that V is 
given by a translation-invariant m-potential. This last condition means that we can 
assume (if necessary by passing to a physically equivalent interaction) that there 
exists at least one configuration u> minimizing all &ox, i-e., 

*ox(a;) = min*ox(w'). ( 2 - 22 ) 

for all X. For any m-potential, the set of all configurations for which Eq. ( |2.22D holds 
is the set of ground states of 3>o- 
(P2) The perturbation interaction T is in some space Banach space B r for some r > 0. 

Since, by condition (PI), the ground states can be determined locally, there is a corre- 
sponding decomposition of the Hilbert space Ha for all A: 

n A = n 1 ^ e n h 2 gh , (2.23) 

where Ti}™ is the subspace spanned by the ground states of V. We can decompose any 
operator Ka € B(Ha) according to their action on Ti 1 ™ and 7i A lgh : 

K A = K A + K\ h + K 1 }, (2.24) 

with 

K A n l x w c n l x w K A n h A igh = , 
Kfn h l gh c n h 2 sh KfnT = , 

K l }n l x w c n\ igh K^n A gh c n l x w . (2.25) 

Accordingly we decompose any interaction T: 

T = T 11 + T hh + T lh (2.26) 

The following theorem shows that, for any integer n > 1, it is possible to construct 
an interaction equivalent to H with the property that is block diagonal up to 



order n. Note that this is a constructive result and a algorithm is given in [DFFR] which 
allows to construct the unitary transformations and the interactions H ^ . 

Theorem 2.1. Consider an interaction of the form 

H = V + XT (2.27) 

where V satisfies Condition (PI) andT satisfies Condition (P2). For any integer n > 1 
there is r n > and A n > such that for |A| < A n there is an interaction = V € 
B Tn , equivalent to H , with 

||T( n ) lh || rn = 0(A n+1 ). (2.28) 

This theorem is useful to analyze the low temperature behavior of quantum spin systems 
when the ground states of V have infinite degeneracy and T lifts this degeneracy (totally 



s 
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or partially) . Consider for example the typical case where the degeneracy is lifted in second 
order perturbation theory. In that case we may take n = 1 and we have T^ 111 = <3(A 2 ): 

flM = V + XTf )n + £ X j Tf )hh + £ XTf )lb . 



(2.29) 



We then decompose 



3>l J 'I J>2 

V + T into a new "classical part" V given by 



2 



V + > X J T 



(i)ii 



(2.30) 



and T contains all remaining terms. The new perturbation satisfy the bounds T = 0(A 3 ), 

T = 0(A), and T = 0(A 2 ). If V is a classical interaction with a sufficiently regular 
zero-temperature phase diagram, then Pirogov-Sinai techniques can be applied to study 
the phase diagrams of V + T for sufficiently small A (see below) . 

Note that this perturbation scheme is not only useful to analyze the low-temperature 
behavior of the model. The new "classical part" V does not need to be classical at all. 
For example, see [DFFR, pFF2j ], if one applies this perturbation scheme to the Hubbard 
model at half-filling, V is given by the Heisenberg model and this gives a rigorous proof 
of the equivalence of both models up to controlled error terms. 

3. Phase diagrams, contour models, and Pirogov-Sinai theory 

A phase diagram in Thermodynamics is a partition of a space of physical parameters in 
domains corresponding to phases; the free energy varies very smoothly inside a domain. 
However, first derivatives or of higher order may have discontinuities when crossing the 
boundary between two domains, and in this case one talks of phase transitions. 



The first proof of a phase transition was proposed by Peierls for the Ising model [Pei]. 
It was extended by Pirogov and Sinai [PS, Sin] to situations where different phases are 
not related by a symmetry. Important extensions and simplifications of the Pirogov-Sinai 
theory include Kotecky and Preiss [KP], Zahradmk [2ah], Bricmont et.al [ BKL| and BS ] , 
Borgs and Imbrie [BI], Borgs and Kotecky [|BK| , BK2 |. An exposition of the Pirogov-Sinai 



theory can be found in EFS]. 



Another extension of the Peierls argument was done in Frohlich and Lieb [FL| using 



reflection positivity ||FSS| , pLSfl . 



3.1. Phase diagrams. We consider the Banach space B r of periodic interactions, with 



the norm defined in (2.6). Here r is any positive number, but further assumptions (bounds 
for the weights of the contours, see below) can be verified in given models only if r is large 
enough. To a given interaction H G B r and temperature (5 we associate the set of all 



translation invariant (or periodic) KMS states or, equivalently [Ara, [si], the set of all 
tangent functionals to the free energy f(H). The set of periodic KMS states forms a 
simplex, so that it is enough to describe the extremal states, or the corresponding tangent 
functionals. We denote the set of extremal states by £@(H). 

In order to define a phase diagram we consider a smooth (p — l)-dimensional manifold 
on the Banach space B r of periodic interactions; it is described by an application u \— * H u , 
from a connected open set U C MP' 1 into B r . For m = 1, 2, 3, . . . , we introduce i?( m ) = 
{H £ B r : \£P(H)\ = m}; accordingly, we partition the set U as 



U 



u u (m) 

m=l 



(3.1) 
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where u E uS m ^ iff H u E E^ m \ The decomposition ( |3.1[ ) is called the phase diagram of 
H" 

The phase diagram of H u , u E IA C M p_1 , is said to satisfy the Gibbs phase rule if the 
following conditions hold. Here, we call "boundary" of «W the set (U® \U®)nU, with 
Z?W the closure of Z/W. 

(i) U = UM U ••• U^ p ). 

(ii) (a) Z^ 1 ) consists of p connected components, each of which is a (p — l)-dimensional 

manifold. The boundary of IA^ is IA^ U • • • U U®. 

(b) ZY^ 2 ^ consists of (!) connected components, each of which is a (p— 2) -dimensional 
manifold. The boundary of is tf( 3 ) U • • • U U®. 

(c) consists of ( q) connected components, each of which is a (p — g)-dimensional 
manifold. The boundary of is Z/« +1 ) U • • • U 

(d) £-/( p ) consists of a single point «o- 

In other words, the phase diagram of H u satisfies the Gibbs phase rule iff it is homeo- 
morphic to a connected, open neighborhood IA' of the boundary of the positive octant of 
MP, in such a way that uq is mapped onto the origin, IA^ P ~^ is mapped onto the union of 
axis Uj{a,, > 0, dj = 0,j ^ i}, and so on... 

Connected components of IA^ are the one-phase region, or pure phase region, IA® is 
the region of coexistence of two phases, . . . , lA^ 1 is the point of coexistence of all p phases. 

We will call a phase diagram which satisfies the Gibbs phase rule regular if the free 
energy is a real analytic function of u in each one-phase region, and if all connected 
components of the manifold IA^' are smooth (C 1 ). 

3.2. Contour models. A contour A is a pair (A, a), where A C 7L V is a finite connected 
set and is the support of A; to describe a, let us introduce the closed unit cell C{x) C W 
centered at x, i.e. C(x) = {y E M? : \y — x[oo ^ i}. The boundary B{A) of A C Z w is the 
union of plaquettes 

B(A) = {C(x) n C{y) :x€ A,y^A}. (3.2) 

The boundary B(A) decomposes into connected components; each connected component 
b is given a label tfy, E {1, . . . , p}, and a = («{,). 

Let AcZ" finite, with periodic boundary conditions. A set of contours {A\, . . . , -4.^} 
is admissible iff 

• Ai C A, and dist (Ai, Aj) ^ ltfi^j. 

• Labels aj are matching in the following sense. Let W = A \ VJ k - =1 Aj] then each 
connected component of W must have same label on its boundaries. 

For j E {1, . . . ,p}, let Wj be the union of all connected components of W with labels j 
on their boundaries. 

For each j E {1, . . . ,p}, we give ourselves a complex function gy u ("free energy of a 
restricted ensemble"), that is real analytic in u E IA. We suppose that the limit (3 — > oo of 
,u exists, and we write 

= lim Reg?> u , 1 < i ^ p, (3.3) 
e£ = mine^. (3.4) 

i 

We consider the partition function ( |2.8| ) for an interaction i3" M = V u + T, where the 
periodic interaction T is a perturbation of V u . We assume that the partition function can 
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be rewritten as 

z a u,t = E U wP ' u,T Mfl e '^' vm ' ( 3 - 5 ) 

{A u ...,A h }3=l i=l 

where the sum is over admissible sets of contours in A.[] The weight w^' u ' T {A) of a 
contour A is a complex function of (3, u, and T, that behaves nicely for (3 large and T in 
a neighborhood of 0. Precisely, we assume that there exists a set W C x B r , that is 
open and connected, and whose closure contains (oo,0); furthermore, we suppose that for 
all u G U and all T) G W, and all contours .A, 

• w@' u,T is periodic with period i.e. we have wA u,T (T a *4) = w^' u,T (A) for all a G 
(£1i) u and all .A. Here r a is the translation operator. 

• |w /3 ' n ' T (.4)| ^ e _/3e ol^l e -1 "'" 4 ' for a large enough constant r (depending on z/, p, and 
I). Furthermore, \^-w^'^ T (A)\ < /3|A|C e~^ e o l A l e - T l A l and \^w^ u ' T+7 ^ K (A)\ < 
/?|^|C||K|| r e _/3e ol A l e _r l A l for a uniform constant C. 

• lim^oo limT^o w^ ,u,T (A) = 0. This means that the weights represent the correction 
to the situation {(3 = oo,T = 0). 

• w^' u,T (A) is real analytic in u; for all K G B r , w /3,u ' T+riK (A) is real analytic in r\ in 
a neighborhood of (the neighborhood depends on K). 

Finally, the free energy is 

fW = -\ lim -LlogZf^. (3.6) 

1 (3 a/Z" |A| A v ; 

We also assume the following properties for fP' u ' T ; 

• fP' u ' T is real, and concave as a function of T; 

• whenever H u + T = JT"' + T', we have 

Although these properties seem difficult to verify in the context of a contour model, they 
are usually clear in the original physical model. 

3.3. The Pirogov-Sinai theory. The results of the Pirogov-Sinai theory are usually 
presented in terms of existence of many Gibbs states for a given interaction. However, it 
is more convenient to think of the Pirogov-Sinai theory as to express the free energy in a 
suitable form for the description of first-order phase transitions: the free energy is given 
as the minimum of C 1 functions ("metastable free energies"), that intersect themselves 
by making angles. Hence a first-order phase transition when varying parameters so as to 
cross an intersection. 

The free energy at zero temperature is given by ( |3.4| ); in typical situations this is the 
minimum over energies of some important configurations (the "potential ground states"). 
The Pirogov-Sinai theory shows that in contour models, this structure extends at low 
temperatures. In the quantum situation one is also interested in adding a perturbation to 
a "nice" model; the metastable free energies then depend not only on /3, but also on the 
quantum perturbation. 

We claim that the Pirogov-Sinai theory allows to construct metastable free energies 
that satisfy the following properties. 



lr The sum includes the case k — 0, and the corresponding term is 53f=i e l3a ' ' A ' • It i s however 



irrelevant, since it does not contribute to the infinite- volume free energy 
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Properties of the metastable free energies. We consider a contour model that sat- 
isfies the structure described in Section 3.&. Then there exist p real functions f^' u,T for 
(f3,T,u) eW xU, such that 

(a) /A«,t = m[Qi f ^ T ; 

(b) L,_li m J;/f*> = e ? , ^^m^W-te; 

(c) for all K £ B r , there exists a neighborhood Mk °f such that fP> u ' T+r i K i s c 1 as 
a function of (u,n) in U x Mk, and |^/f ,u,T+vK \ <C (7||i^|| r for a constant C 
depending on u, p, £ only; 

(d) fi' u ' T is a real analytic function of u in VJl^y = {u : /f'"' T < fj' u,T Vj ^ i}. 

Notice that the point (d) implies that the free energy fP' u > T is a real analytic function 
of u in Uj97l{j} (which is the region of uniqueness, as will be seen below). 

The proof of these properties involves the full artillery of the Pirogov-Sinai theory. The 
item (c) is not really standard and may appear as superfluous technicalities, but it plays a 



role when establishing the properties of the phase diagram, see Theorem below. Since 
the present paper is only aimed at studying a special class of quantum models, we content 
ourselves with an outline of the proof, so as to make it plausible for readers who have 
knowledge of the details of the Pirogov-Sinai theory. A review of the Pirogov-Sinai theory 
is expected to appear shortly and will contain a detailed proof of these properties. 



Sketch of the proof of these properties. We heavily rely on [ BKU |, which itself follows [PS 



Sip , Zah, BI, BK, BK2f |. Our metastable free energies are defined as the real part of the 



metastable free energies of [BKU], which are complex in general. 

The first step consists in defining the metastable free energies. This can be done by 
introducing truncated contour activities and truncated partition functions following the 

(n) 



inductive procedure of flBKU ], Eqs (5.6)-(5.12). One obtains metastable free energies /, 



j 



(that depend on f3, u, T). One can then prove the claims of Lemma A.l i), iii), iv), v), vi) 



of iBKUj . We set then /f' n ' T = li m7W0O & ] . 



At this point we have well-defined metastable free energies depending on /3, u and T 
(that is, they are functionals on the Banach space of interactions), and the free energy of 
the system is given by the minimum of the metastable free energies, as stated in item (a). 
It is also clear that linig^oo limi^o fi' u,T = ef, and that /■ U ' T is real analytic in u on 
5Dt{j}. What remains to be done is to check differentiable properties. 

For given T and K, we consider fP' u ' T+1 i K a s a function of {u,n). This is a mild 



complication of the situation in | BKU |, since the metastable free energies here depend 
on p parameters instead of p — 1. One then gets the items ii) and vii) of Lemma A.l 
- the partial derivatives with respect to 77 of the truncated contour activities and of the 
partition function with given external label satisfying the claims of the lemma with a 
constant Co||-?f|| r instead of Co- 

Finally, the metastable free energies are given as convergent series of clusters of contours, 
the weights of those obeying suitable bounds. This leads to item (c). □ 

We show now that these metastable free energies allow for a complete characterization 
of tangent functionals, under the extra assumption that the situation at zero temperature 
and without perturbation satisfies the Gibbs phase rule in a strong sense. 

The stronger condition for the Gibbs phase rule is that, for some uq 6 U, we have that 
all "potential ground state energies" are equal, e"° = e"° for all i, j, and that the matrix 
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of derivatives 

(w-H- e p]) ( 3 - 8 ) 

\OUj 1 PJ /1 ^p-l 

has an inverse that is uniformly bounded. Actually, energies e" may not be differentiable; 
in this case, we consider the same matrix with Regf' M instead of ef, and we suppose that 
it has an inverse for all (3 large enough, the inverse matrix being uniformly bounded with 
respect to u 6 U, and (3 ^ const. 

Theorem 3.1 (Stability of the phase diagram). Assume that there exist metastable free 
energies fi' T , 1 ^ i ^ p, that satisfy all points (a)-(d) of the properties above. We 
assume in addition that the strong version of the Gibbs phase rule, described above, is 
satisfied. 

Then for (3 large enough and \\T\\ r small enough (depending on p and on the bound of 



the inverse of the matrix of derivatives (3.8) ), there exists U 1 C U such that the phase 



diagram for H u + T , u E W , at inverse temperature [3, satisfies the Gibbs phase rule and 
is regular. 



Theorem states that there exists u' Q € W such that the set of tangent functionals 
to the free energy at H u o + T is a simplex with p extremal points. More generally, we 
have the decomposition W = U'^ 1 ' U • • • U U'^ such that for u £ W' , the set of tangent 
functionals at H u + T is a g-dimensional simplex. 



This "completeness" of the phase diagram was addressed in [Zab] and [BW]. The 
approach was however different and involved studying the Gibbs states, which is more 
intricate and does not easily extend to the quantum case. It is simpler to look at tangent 
functionals, and then to use existing results on their equivalence with DLR or KMS states. 

Notice that the Pirogov-Sinai theory also provides various extra informations, such as 
the fact that the limit of W , as T — > and (3 — > oo, is equal to U^ q \ Also, the 
extremal equilibrium states can be shown to be exponentially clustering. We do not claim 
these properties here however, because doing so would require extra assumptions and 
technicalities in the description of the abstract contour model. 



Proof of Theorem 3.1. Items (b) and (c) of the properties of metastable free energies (with 
rj = 0) imply that there exists u' such that f^' u °' T = j^ u o< T f or a rj ^ ^ anc [ th e ma trix 
of derivatives 

has a bounded inverse, uniformly in it in a neighborhood W of u' . Let us define 

Wli = {ueW: t u ' T = min J/'"- 7 '}. (3.10) 

and, for Q C {1, . . . ,p}, 

m Q = p| m, t \ |J wti (3.ii) 

(notice that SPt/j) C SJtj). By the implicit function theorem, each is described by a 
C 1 function from an open subset of RP - '^ into U' . If we set = U|q| =(J 9JTq the phase 
diagram satisfies the Gibbs phase rule, provided there are exactly \Q\ tangent functionals 
at H u + T for each u£Wt Q . 



QUANTUM LATTICE MODELS AT INTERMEDIATE TEMPERATURE 



13 



Each metastable free energy fj' u,T , j E Q, defines a tangent functional ay: for all 
K E B r , we set otj{K) = ■^fj ,u ' T+r,K \n=o- Notice that item (c) ensures boundedness 
of the tangent functional.^ We show now that these tangent functionals are linearly 
independent, and that any other tangent functional is a linear combination of these ones. 

We examine the manifold where q phases coexist; without loss of generality, we can 
choose u E 9Jtg with Q = {1, . . . ,q}. The determinant of ( ^S| ) can be written as a linear 
combination of determinants of 

d .[f^-f^) (3 . 12) 



with k±, . . . , fcg-i being q — 1 different indices. Since the determinant of ( |3.9| ) differs from 
0, at least one of the determinants in the previous equation differs from 0. Without loss 
of generality we can assume that 

is not singular. 

Our analysis is local, so we can take u = and H u = H° + YJj=\ u j K j- Then 



(3.7) implies that otj(Ki) = -g^-fj ,u ' T \ u =o, and non-singularity of ( p.!3[ ) shows that 



a. 



■]■> 

1 ^ 3 ^ <Z> are linearly independent. Furthermore, it also implies that for all tangent 
functional a 1 the system of equations for £ = . . . , £ q ), 

q 

a/(K i )=J2$ j a j (K i ), i = l,...,q-l, (3.14) 

j=l 

has a unique solution with ^ • £j = 1 . Now we consider any K £ B r ; we define gj(u, rj) = 
f^ T+ " K , 1 < j < «, and 

(5l(«,»7) - 0g \ 
: . (3.15) 

9 q -i{u,rj) -g q (u,rj)J 

We have g(0, 0) = 0, Jj:<?(0, 0) is an isomorphism, and g(u, rj) is a map of class C 1 by item 
(c) of the properties of page |ll[ By the implicit function theorem there exists a map u(rj) 
such that g(u{rj),r]) = 0. We introduce the interactions 



1 9-1 

R( V ) = K+-^2u j (f,)K j , (3.16) 



<z-i 



= lim #(77) = K + > n' (0)Kj. (3.17) 

J'=l 

Then using ( |3.7| ) we have 

f/3,0,T+riR(n) _ fPfi,T+r)R(r]) _ _ y/3,0,T+»jR(r;) ^ 



2 One may wonder whether the functional ay is linear. It is actually, because Qj can be obtained as the 
limit of linear functionals that are tangent to the free energy, uniquely defined for all points of ?fft{j} — a 
region of parameters where the concave free energy has a unique tangent functional. 
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Differentiating with respect to rj, we obtain (recall that a' is tangent to f/ 3 '°> T+r iR( , n) a t 



a'(R) = a t {R) 



a q (R). 



(3.19) 



Then obviously a'(R) = £^ ■ £ja.j(R), and it follows by linearity of the tangent functionals 
that 



3=1 



(3.20) 
□ 



4. Results of the quantum Pirogov-Sinai theory 



We summarize in this section the results obtained in |BKU, DFF, DFFR, KU], and in 
the present paper. All results concern the situation where the interaction has the form 
H = V + T, where V is a classical interaction satisfying the standard Pirogov-Sinai 
framework, and T is a small perturbation. The temperature will be assumed to be small. 
The results however split into four classes, according to whether we use the perturbation 
methods of [ DFFR[ ] (Section 2.2), and whether we include high temperature expansions 
to analyze phases at intermediate temperatures. 

In this section, we implicitely assume all properties of the metastable free energies, see 
page |ll|, to be valid — without these properties the statements below would not include 
completeness, i.e. we could not ascertain to have identified all the periodic Gibbs states of 
the systems. 



4.1. Quantum perturbation of classical model with finitely many ground states. 

In this case the classical interaction V has finitely many ground states and the phase 
diagram of V + T is, at low temperatures and for sufficiently small T a small deformation 
of the zero temperature phase diagram of V . The extension of the Pirogov-Sinai theory 



to this class of quantum systems goes back to [ Pir | and was proved in [BKU, DFF]. 



(a) Structure: We denote by fi = {1, . . . , M} the space of classical configurations; the 
dimension v of the physical space is always supposed to be bigger or equal to 2. The 
interaction has the form H = V + T, where V is a block interaction and is diagonal with 
respect to the basis of classical configurations: if A = U(x) = {y : \y — x\oo ^ R} for 
some x G Z y , 

V A \e) u = $x(uu(x)) (4.1) 

and Va = if there is no x with U(x) = A. The function & x depends on fi 6 W C M p_1 , 
and we assume that its derivatives ^f^x(^u(x)) are bounded uniformly in x,/u,,u;,j. 

A finite set G = {g^\ . . . ,g^} C of periodic configurations is given, that contains 
all ground states of V for all fi (see below the precise assumption). We write Ga = {oa '■ 
g G G}. We suppose that $ x {du(x)) 1S independent of x, for all g G G, and this value is 
denoted by eg (this is the mean energy of the configuration g). 
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(b) Assumptions: 

(Al) A gap separates the excitations: for all usjj(x) £ Gu(x)> 

(uniformly in fj,). 

(A2) The zero temperature phase diagram is (linearly) regular: there is n € U such that 
eg = e^,° for all g, g' & G, and the inverse of the matrix of derivatives Mq, see ( |3.8|) , 
is uniformly bounded. 

(c) Properties of Gibbs states: 

Theorem 4.1. Assume (Al) and (A2) hold true. There exist Po,c < oo (depending 
on v,R,p,M and on the periods of {g^'} and H only) such that if (3 A ^ /?o and 
||T|| C /A ^ 1, the phase diagram of the quantum model satisfies the Gibbs phase rule and 
is regular in a neighborhood hi' Cl/ of n . 

In the single phase region, i. e. if /j£ m^({g}), the KMS state w^ T {-) is close to the 
ground state g: for all K € Oa, nm /3->oo,||T|| r ->o w^'^' T {K) = (e g \K\e g ). 

The condition ||T|| C /A ^ 1 means that T is a perturbation with respect to V\ c plays 



the role of the perturbative parameter: from the definition (2.6) of the norm || • || c , \\Ta\\ 
must be very small if c is very large. 



The proof of this theorem follows from [ BKU , DFF ] . 



4.2. Models with infinite degeneracy. Consider a model whose classical part has 
infinitely many ground states, and a perturbation which lifts this degeneracy completely. 



The pertubation methods of [DFFR|(see section 2.2) permits in certain cases to analyze 
this by constructing an equivalent interaction with a new classical part which has finitely 
many ground states. In this case the new perturbation has a slightly more complicated 
form than in Section 4.1 and the following theorem deals with this situation. This situation 
was considered in [ DFFRj (for a different approach see | KUf ). 



(a) Structure: The space of classical configurations is again = {1, . . . ,M} Z ". We con- 
sider two sets G, D C n, with D C G finite, D = {d^, . . . , d^} is a finite set of periodic 
configurations; G may be infinite and will represent the configurations of low energy. For 
Ac2", the Hilbert space Ha has the following decomposition Ha = H l ™ © H A 1&h where 
H 1 ™ is the subspace spanned by the low energy configurations gA € Ga- The interaction 
has the form H = V + T, where V is a classical block interaction with uniformly bounded 
derivatives -gjj--<& x (uju^), and T is a perturbation that is submitted to some restrictions, 
see the assumptions below. 

(b) Assumptions: 

(Bl) A gap separates high and low energies: for all uj\j{x) ^ Gjj^, 

®x(vu(x)) -™g$x(9u(x)) > A o- 

(B2) Gap with the ground states: we assume that $x(djj(x)) 1S independent of x for d G D, 
and for all u V ( x ) £ D u( ^ x) , 

®x(wu(x)) -*P&®x(du(x)) > A 

(and we assume that A ^ Ao). 
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(B3) The perturbation may be decomposed T = K + K' + K"; for all A, 
K' A U l T = 0, K' A H h 2 gh cH* igh - 



A 



(there is no assumption on K)^ 
(B4) The zero temperature phase diagram is (linearly) regular, i.e. all energies are equal 
for some /x G and the matrix [see ( |3.8|) 1 has a uniformly bounded inverse. 

^ Properties of Gibbs states: 

Theorem 4.2. Assume (B1)-(B4) hold true. There exist Po,c < oo (depending on 
u,R,p,M and on the periods of{d^} and H only) such that if (3 A ^ (3q, \\K\\ c /A ?J 1, 
||JiC || c /Ao ^ 1, \\K"\\ C / Aq sj 1 the phase diagram of the quantum model satisfies the 
Gibbs phase rule and is regular in U' C U, W 3 /i . 

In the single phase region, i.e. iffie Wt^dd}), the KMS state w p ^' T {-) is close to the 
ground state d: for all K £ Oa, li m /3-*oo,||T|| r -^o w^ ,fJ/,T (K) = (ed\K\ed)- 

The proof of this theorem is given in [DFFRj. A somewhat different method yielding 



similar results has been developed later in [ |KU| |. 

4.3. Combined high and low temperature expansions. Here we consider models 
whose classical part V has partially ordered ground states, typically described by peri- 
odic configurations of holes and particles but still with infinite degeneracy due to, e.g., 
degeneracy of the spin at each site. Together with the quantum perturbation the system 
may have a continuous symmetry. We will suppose that the temperature is low and, in 
addition, that /3||T|| C is actually small (i.e. the temperature is large compared to T) and 
we will prove that in this case one phase corresponds to each periodic configuration of 
holes and particles and that in this phase the spin degrees of freedom are in a disordered 
phase. This situation has many similarities with that of [BKL|, and could be called "a 



theory of restricted ensembles in quantum lattice systems". 

(a) Structure: As before, let £1 = {1, . . . ,M} Z ". Intermediate temperature phases will be 
characterized by "motives" giving partial information on the underlying configurations. 
In order to describe this, we consider a partition of {1, . . . , M}: 

JV 

{1,...,M} = (J Ij with k nlj = $. (4.2) 

We denote A/a = {1, . . . , N} A (and A" = N%»). For n G A", we write Q n = {u G O, : u x G 
I Ux Vx}. Let G = {g^ l \ . . . ,g^} C M be a finite set of periodic configurations; this is the 
set of motives and a pure phase will be associated with each of these configurations. We 
write VLq = U ge c^g- 

The interaction has the form H = V + T, where V is a classical block interaction with 
uniformly bounded derivatives w.r.t. and T is a perturbation. We introduce restricted 
partition functions for each g G G: let 

Z 9 K = e - ^^^)^**^ 3 ^ (4.3) 



3 Motivation comes fro m (2.25). It is however slightly more general, and it is just what is required in 
the proof of Theorem 4.2. 
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and 



h^ = -\ lim r^logZl (4.4) 

9 P k/ir |A| A 



The ground energies are eg = linig^oo hg ,,J/ , g G G. 

(b) Assumptions: 
(CI) For all configurations ujjj^ ^ ^g,U{x)i we have 

®x{uu(x)) ~ min ®x{u' u{x) ) > A. 

Moreover, we assume that 

min <f>x(uu(x)) = e M (ff) 

independently of x, for all g G G. 
(C2) We need a condition that ensures that no phase transition takes place in a restricted 
ensemble fl g ; in other words, spatial correlations should decay quickly enough. The 
following condition is stronger, and amounts to saying that there is no correlation 
between different sites. For all g G G, we suppose that there exists an on-site 
interaction <I> 9 such that for all x: 



for all u) G £l g . 



(C3) The zero temperature phase diagram is regular with eg = e^P, g,g' G G, for some 
fj,Q G U, and the matrix Mq, see (3.8), has a uniformly bounded inverse. Q 



(c) Gibbs states at intermediate temperature: 

Theorem 4.3. Assume (C1)-(C3) hold true. There exist (3q,c < oo (depending on 
u,R,p,M and on the periods of {g^} and H only) such that if Pq ^ /3A < oo and 
/3\\T\\ C ^ 1, the phase diagram satisfies the Gibbs phase rule and is regular in hi' C IA, 

In the single phase region, i. e. if H G m^dg}), the KMS state w (3 ^ T (-) is close to the 
motive g: for all K G O a , lim /9 _ +00) || T ||,._,o wP>»> T {K) = (Tr(P A ))- 1 Tr(ErP A ), where P A 
is the projection given by E WA en A \ e "A)( e - 



Remark: It follows from our assumptions that T is small compared to V; more precisely, 
||T|| C /A < 1/Pq. 



This theorem is actually a consequence of Theorem 4.4 below, see the remark on page 



19 



4.4. Infinite degeneracy, high and low temperature expansions. Here we consider 
systems where phases result from subtle interplay between potential and kinetic energy, 
combining the effect described in Section 4.2 and 4.3. The quantum perturbation lifts par- 
tially the degeneracy of the classical interaction, leading at intermediate temperatures, to 
spatially ordered phases. Hereafter we describe the general framework in a rather abstract 
way; it will be illustrated in Section [|, and the reader may gain better understanding by 
working out a concrete application. 



4 If {eg } are not C 1 , we consider the matrix of derivatives of hg'^ for f3 large; it must have an inverse 
that is bounded uniformly w.r.t. /i and large f3. 
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(a) Structure: The space of classical configurations is fi = {1,... ,M} Z "; we consider a 
partition like in ( |4.2| ) and define similarly M and £l n . We consider a (possibly infinite) set 
G C N that represents low energy configurations; the Hilbert spaces decompose in the 
following way: Ha = H l ™ ® T^-a^ > wnere H l ™ is the subspace spanned by the low-energy 
configurations qa £ Ga- The interaction has the form H = V+T; V is a block interaction 
with uniformly bounded derivatives ■$^:$x{uu(x))i the perturbation T decomposes further 

T = K + K' + l^"; we shall require different assumptions on K , K' , K" , motivated by 
the perturbation theory of Section 2.2. 

We suppose that a finite set D = {rfW, . . . , d^} C G is given, that corresponds to 
possible ground states. For each d £ D, we define the corresponding restricted partition 
function 

Z d k = ^2 e^^^O-OcA**^*)) (4.5) 
and the corresponding restricted free energy 

and = lim^oo h^. 

(b) Assumptions: 

(Dl) A gap separates high and low energies: for all ^u(x) & ^G,U{x)i 

$x(uJu(x)) ~ max $ x (u' u(x) ) > A . 

(D2) Gap with the ground states: for all oju( x ) £ &d,u(x)i 

(D3) For all d £ D, there exists an on-site interaction $ rf such that for all u € 0,$ and all 

x, 

^(%( I )) = ^K). 



Moreover, we suppose that 



min (a;,.) = e£ 



independently of x. 

(D4) The quantum perturbation T = K + K' + has same properties as in (B3), with 
respect to the decomposition into low and high energy states. 

(D5) There is /x £ U such that = effl, d,d' € -D, and the matrix of derivatives ( |3.8[ ) 
has a uniformly bounded inverse (see the footnote of (C3) if is not differentiable). 

(c) Properties of Gibbs states: 

Theorem 4.4. Assume (D1)-(D5) hold true. There exist (3q,c < oo (depending on 
v, R,p, M and on the periods of {d^'} and H only) such that if (3q /3A < oo, /3||1C|| C ?J 
1, H-f^'llc/Ao ^ 1, ||-K""|| c /Ao ^ 1, and f3\\K"\\^/ Aq ^ 1, the phase diagram satisfies 
the Gibbs phase rule and is regular in an open setW C IA that contains /i . 

In the single phase region, i.e. if fie WiP({d}), the KMS state w^' T {-) is close to the 
motive d: for all K G O a , ^-p-^< x ,\\T\\r^ w ^ ,T i K ) = (Tr^)) -1 ^^-^), where P A 
is the projection given by T,cj A en d , A \^a){^a\ ■ 
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This theorem follows from the contour representation obtained in Section ^, together 
with the Pirogov-Sinai theory. 



Remarks: 1. Theorem 4.3 is an immediate consequence of Theorem 4.4. Indeed, we 
clearly recover the setting of Section |4,3j by choosing G = O (i.e. all configurations have 
low energy), and K' = K" = 0. 

2. These two theorems also generalize results of [ Uel | : they can be applied to the 
Hubbard model 



H 



-t Y (4<7<V + h.c.) + U J^n^TVi, 

<x,y> x 



(4.7) 



to show that the high temperature phase extends to 

{ (p, t, U) : pt small} and { (fi, t, U) : (3t 2 /U small} 
(standard high temperature expansions apply when both fit and f3U are small). 



5. Example: Extended Hubbard model 

This is a Hubbard model where particles interact among each other when their distance 
is smaller or equal to 1. Explicitly, 

( C L<V + h - c + U Y n ^ n *l + W Y UxT 

xGA <r,j/>cA 



-t 



£ 

<x,y>CA 



nY n * 

xGA 



(5.1) 



Here, ct,a, c xa are creation, annihilation, operators of a fermion of spin a at site x; <x,y> 
stands for a set of nearest neighbor sites; n xu = c xa c xa is the number of particles of spin 
a at x (it has eigenvalues and 1); n x = n x ^ + n x \^ is the total number of particles at x. 
The coefficient t represents the hopping, and will be taken to be small compared to the 
nearest-neighbor repulsion W; /u is the chemical potential. The classical limit t — > was 
studied in [ Jed , BJK|. The stability of the chessboard phase 50t(o,2) ( see below) with small 
t is a straightforward application of [DFF]; a later study devoted to it is [ |BK3[ . 

We start by analyzing the classical interactions. The configuration space is Q = {0, j, J. 
, 2} z " and the corresponding classical interaction can be written as (taking R = ^) 

U 



y<EU(x) 

Here we introduced q y € {0, 1, 2}: 



x - W 

2^ d ^. 2 + ^T 



<y,z>cU(x) 



2 ;/ 



y£U(x) 



(5.2) 



if LVy 
if OJy 

if uj v 



-- 

-] or uj y 
-- 2. 



(5.3) 



The interaction can also be written as a sum over pairs of n.n. sites; this simplifies 
the analysis of the zero temperature phase diagram, and the search for symmetries (see 
below). This pair interaction is given by 

U 



®<x,y>{qx, Qy) = ^(<W + 6q y> 2) + Wq x q y - -^-{q x + q y ). 



(5.4) 



This model has a hole-particle symmetry. Introducing the unitary operator U such that 



ucLu- 1 



and U c xa U 



c xa , we see that UTfJJ 1 = T\. As for the potential, the 
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effect of the symmetry can be exhibited by considering classical configurations; defining 
q' x = 2 — q x , and fjf = U + ^vW — fj,, we easily check that 

= ty) + C ( 5 - 5 ) 

where C = —U /u—AW+2^/u does not depend on (q x , q y ). As a result, the phase diagrams 
(U, ji) are symmetric along the line 



(5.6) 



for any temperature. 

The zero temperature phase diagrams with t = are depicted in Fig. in both cases 
< and W > 0. 





u\W\ 


-2 




M 



u 




(a) (b) 
Figure 2. Zero temperature phase diagrams of the extended Hubbard 
model, (a) when W < and (b) when W > 0. The dashed line represents 
the hole-particle symmetry, see (|5.6|). 



In the case W < 0, it decomposes into three domains TIq, SOTi, and 97^2 ; 9^0 and 9JT2 
have a unique translation invariant ground state with respectively and 2 particles at 
each site. In SDTi, any configurations with one particle per site is a ground state; there is 
degeneracy 2' A ' since each particle has spin f or j. 

The situation W > presents a richer structure with six domains. Domains 9Jto> 
and 5D?2 have same features as with attractive n.n. interactions. In between appear now 
domains SPt(o,2) > ^(1,2) an d 9^(o,i) • ^(0,2) consists in two ground states, the two chessboard 
configurations with alternatively and 2 electrons per site. 50t(o,i) has 2 • 2al A l ground 
states of the chessboard type, one sublattice being empty, while the other has exactly one 
particle of spin f or j; 971(1,2) is similar, with 2 particles per site on one sublattice and one 
on the other. 

We are interested in the case where the temperature is small, but bigger than 0, and 
with small hopping. The phase diagrams for large (3 and small (3t are presented in Fig. ^. 

In the case W < 0, all three domains survive at low temperature and with t / 0; a 
first-order phase transition occurs when crossing the border between any two domains. 
The point (-Mr? = 2, = 1) belongs to 97t^'*: this phase has residual entropy (it also 
has more quantum fluctuations, although this has much less effect). The Gibbs state 
corresponding to the domain 971^'* is thermodynamically stable and exponentially cluster- 
ing. The restriction to intermediate temperatures (f3t ^ e) is important, because, for 
v ^ 3, a phase transition is expected when the temperature decreases, leading to an 
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v\W\ 

2 / 




/ 4 


-2 


<•* 



SOT*' 



u 

u\W\ 



f3,t 



U 



(3,t 



0Kg 

(a) (b) 
Figure 3. Phase diagrams of the extended Hubbard model at intermediate 
temperature and with small hopping, (a) when W < and (b) when W > 0. 
Bold lines denote first-order phase transitions. White is the region V e that 
resists rigorous investigations, where second-order transitions are expected. 

antiferromagnetic phase, that breaks both symmetries of translations and of rotations of 
the spins. 

The phase diagram at finite (5 and non zero t is especially interesting for W > 0. There 
are not six, but only four domains SPt^'*, £0tf*, QJlf'* and Wl%£; see Fig. g. Indeed, the 
three domains corresponding to chessboard phases have merged into a single domain (this 
was first understood and proven in [BJKj] in absence of hopping). The free energy is real 
analytic in the whole domain £0Cu* • The transition between £0?2 and £0tfb is presumably 
second-order, but our results do not cover the intermediate region between these domains. 
The boundary between 93^* and 971^'* contains a part where a first-order phase transition 
occurs, that can be rigorously described. Crossing the boundary elsewhere presumably 
results in a second-order transition. Due to the thermal fluctuations, the segment from 
(2,2) to (2,4) belongs to SUlf *. 

Our results for this model are summarized in the next two theorems. 

Theorem 5.1 (Hubbard model with attractive n.n. interactions). Let v ^ 2. There 
exist constants [3q < oo and Eq > (depending on v) such that the phase diagram 
(U,fi) for f3\W\ ^ /?o and (it ^ Eq is regular; domains 9Jla'*, a € {0,1,2} satisfy 



lim 



linn,- 



9Jt a . If (U, n) belongs to a unique Wlfr*, there is a unique Gibbs 



state. Furthermore, the density of the system is close to a, 

\{n x ) -a\ ^ e{fi,t), 
for all x. e(P,t) can be made arbitrarily small by taking (3 large and t small. 

In order to describe the situation W > we first introduce the region of the phase 
diagram V e where we have no results. Let 



C 



(£0l ( o,2) U 9Jt(i, 2) U 971(0,1) 



f|(aJToU9JTiU£0t 2 )] \ [£0t ( o,2) n 9JTi] , (5.7) 

(5.8) 



and for e > 0, 

V e = (J B e (U,fi) 

where B e (U, fj) is the open ball of radius e centered on (U, fx). We restrict our considerations 
to the complement of V e . 
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Theorem 5.2 (Hubbard model with n.n. repulsions). Let v 2 and e > 0. There exist 
constants (5q < oo and eq > ( depending on v and e ) such that if (3q j3W < oo and 
[3t ^ £q, we have the decomposition 



v c e = sutg'* u m 



U OTtf * U 9Jt^ b , 



and 



are domains with a unique Gibbs state. 



Densities are close to 0, 2, 1 respectively in the sense 



(n x ) ^ e(J3,t) 
(n x ) > 2-s(p,t) 
\(n x )-l\ < e(j3,t) 



in Wl\ 
in 971 
in Wl 



o 

~(3,t 

2 

lit 



(5.9) 



with e(P, t) arbitrarily close to if (3 is large and t small. 

(ii) 9Jtfb* C 971(0,2) ^ ^(1)2) U ^(o,i) ^ s a domain with two extremal Gibbs states of the 
chessboard type. The free energy is a real analytic function of j3 and \i in the domain 

: p /W < < eo/t and (U,y) E 



(in) an&f n is a line of first- order phase transition, with exactly two three extremal 
states. 



7 /3,t 



Remarks: The proofs of Theorems 5.1 and 5.2 use Theorem [Oj. But using Theorem 4.1, 
one could establish stability of domains 97to> a7t 2 > 971(0,2) for all ft\W\ ^ fio, without the 
restriction that the temperature be not too small. Another possible improvement, for 



U, W > 0, would use Theorem iA to replace the condition /3t ^ Eq by (3t 2 /U ^ £q. The 
later clearly allows lower temperatures.^ 



6. Combined high-low temperature expansions 



In this section we simultaneously perform a low and a high temperature expansion. The 
temperature is low, in such a way that excitations above the low energy states (Tl}™ ) are 
rare. At the same time, the temperature is high relatively to the quantum perturbations 
K and K". These expansions allow to write the partition functions as one of a contour 
model, that can be treated by the Pirogov-Sinai theory, see Section 3.2. 



We rewrite the quantum model as a contour model, by making a mixed low and high 
temperature expansion (Section |6.1| ); we define suitable weights, so that the partition 
function takes the form required in Section 3.2. Section ET2] is devoted to proving that 



the weights are small compared to their size. Finally, we explain in Section |6.3j how other 
requirements of Section are fulfilled. 



6.1. Expansion of the partition function. Our intention is to expand in K + K' + K"; 
in order to simplify the notation, we introduce B = {B,i), B C i = 1,2,3, and we 
write K B = T B with B = (B,l), K' B = T B with B = (B,2), and K' B = T B with 
B = {B,2>). We refer to B as a transition. 



5 Furthermore, the restriction to intermediate temperatures arises because of possible antiferromag- 
netism due to "quantum fluctuations" of strength t 2 /U; it should be stable for f3t 2 /U > const; therefore 
this new condition is qualitatively correct. 
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Using Duhamel's formula, we obtain 

Tr e"^ A = Tr e~^cA Vb + ^ £ £ /" dn . . . dr m 

m > 1 Bi,...,B m a)^,...,w™ ■™<TL<...<T m </J 

At this point, it is natural to define the supports of contours as all sites that belong to 
UjBj, or for which there exists cj 3 such that ^ ^D,u(x)- But two technical difficulties 

arise: are periodic rather than translation invariant; and the weight of a 

contour should not depend on the configuration outside of its support (but it may depend 
on the labeling a). The later difficulty is specific to systems with phases given by a 
restricted ensemble instead of a single configuration. To account for these difficulties, we 
introduce a partition of the lattice into cubes of size £, where I is the 1cm of the periods 
of {d^} (considering all spatial directions). 
Let B = \J x& bU{x); we define excited cubes. 

• A cube C is quantum excited if there is Bi such that C D Bi / 0. 

• Otherwise, it is classically excited if there is lo 3 and x € C such that ^ ^D,u{x)- 

Consider the set Q of quantum excited cubes, the set £ of classically excited cubes, and 
the set J\f of cubes that are neighbors of Q U £ (two cubes C ^ C are neighbors iff there 
exist x £ C and y £ C 1 with \x — y\oo = 1). Connected components of Q U £ U N form 
the supports of the contours. Connected components of the complement of Q U £ U M are 
characterized by a configuration d € D, and this information may be stored in the labeling 
a. The union of all components corresponding to d is denoted Wd- Then 

A= Qu£U Mu{U deD W d ), (6.2) 

see Fig. |3| for illustration. Wd is a union of cubes, each cube C contributing in ([O]) by a 
factor [we use (D3)] 

£ e-/9E. 60 *S(<-) = e -K"^\ (6.3) 
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Figure 4. The space A is divided into cubes; contours are formed by 
excited cubes (in black) and by their neighbors. There are four contours 
on this picture. 
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Summing first over admissible sets of contours {Ai, . . . ,Ak}, we can rewrite ( |6.1| ) in 
the following way, 

Tre -/3H A= [JJ e _/3E ^^^ K) 

{Ai,...,A k } deD u>w d £ttd,w d 

[[EE E E / dn...dr m 

&iCQ j j 

e -n E-ex, ( ^ ! ^ ^ } (cj ™ ! ^ ^ } e -0>-r») £^ M*£ w ) (6 4) 

We used here the fact that the contribution of different contours factorizes. There are 
several restrictions to the sums over transitions {B{\ and configurations {a>^ }: each cube 

of Q is intersected by at least one Bi ; {u A . } are compatible with the labeling ay ; excited 
cubes of Aj \ Q do not touch the boundary of Ay, and non excited cubes in Aj have 
at least one neighbor that is excited. In the last line appears configuration with 
U{x) n Wd / 0, hence depending on oJw d - However, in such belongs to a cube 

that is not excited, so that ^i/tx) ^ ^d,u(x)- From (D3) we can substitute ^x(^um) with 

^(u>i), which does not depend any more on the configuration outside the support of the 
contour .0 Then we obtain 

k 

TV e"^ A = E II e-^S"*|Wk| f[ z ( Aj ) (6.5) 
{Ai,...,A k }d&D j=l 

where the sum is over admissible sets of contours, and z(A) is the weight of the contour A. 
The explicit expression for z(A) looks rather tedious, but the main point is to establish 
the properties of Section |3.2[ The expression of z(A) is 



*W= E E E E / dn...dr, 



4 ' 

B.CQ 



e 



-n E,eA*«K(«)) ( u \\T Bl \u 2 A )... (uj%\T Bm \oj\) e'^'^ *«(^ (b) ) j (66) 
with some restrictions on the sums over and {w^}, see above. 

Remark: we constructed contours out of cubes, while the supports of contours in Section 



3.2 are any connected sets. There is no contradiction, if we define the weight z(A) to be 



if the support of A is not a union of cubes. 

6.2. Bounds for the weights of the contours. We turn to the proof of the exponential 



decay of the weight of contours, as required in Section [Sj. We give the following "space- 
time" interpretation to the collection of sums and integrals in ( |6.6D : we view (Bj,Tj) as 
a subset of A x [0,/3] per , with periodic boundary conditions along the "vertical" interval 
[0, 0\. Furthermore, to each "time" r £ [0, /3] per corresponds the configuration lo 3 for 



This is why cubes that are neighbors of excited cubes need to be considered as part of contours. 
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which (rj-i,Tj] 3 t. We define 

m 

B=(jB j x{T j }\jQx{0}; 



3=1 
m+1 



m+1 



E=\J E{J Q ) x [73-1,73], |£| = J] \E{<J Q )\{Tj ~ Tj-i), 
3=1 3=1 
(with ro = 0, r m+ i = (3, and a/™ +1 = oj 1 ). Here, we set 5 = U x& bU(x), and E{ujq) 

{x € Q : ^ f]G,Cf(a;)}' 

/,From assumptions (Dl) and (D2) we can bound 



m > QcA Bi B m 0,1 170 



dri . . . dr„ 



0<Tl<...<T m </3 



-i9A|£|/<- 



e-^IJll^H, (6-7) 



3=1 

where the sums over {Bj} and {uj^} satisfy the restrictions explained above. 

We view each Bj as a connected subset of W +1 (one can e.g. add links between nearest 
neighbors). Then B U E is a subset of 1K^ +1 made out of vertical segments and horizontal 
sets. We consider connected components of B U E. For a connected component with m 
horizontal sets and m! m — 1 vertical segments, we deleted m! — m + 1 of the latters, 
in such a way that the component remains connected. One of these components contains 
Q x {0}, possibly with extra vertical segments and horizontal sets. Other components 
have m horizontal sets and m — 1 vertical segments. Because of the structure (D4), or 
(B3), components not linked with Q x {0}, either consists in a single transition of type 
K, or include at least two transitions of type K or K" . 



P 




A 

Figure 5. A "gather" with 6 transitions and 5 vertical segments. 

A connected object with m horizontal sets, and (m — 1) vertical segments that end on 
the horizontal sets, is called a gather and is denoted by the letter Q. It is illustrated in 
Fig. [|. We introduce the following sets of gathers: 

• <3 m : gathers with m horizontal sets, one containing the origin {x = 0} x {r = 0}; 

= 000. 
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• gathers of ©i that consist in a unique transition of type K. 

• &m- gathers of (S m , with at least two transitions of type K or K"; <& = 0^,. 

The connected component of B U E that contains Q x {0} can be viewed as a set of 
gathers, each gather being connected to Q x {0} by a vertical segment. 

Since a choice of {B{\ and {wg} leads to a set of gathers, we obtain a bound by first 
integrating over sets of gathers, then summing over compatible space-time configurations 
ua, and choosing which gathers are linked to Q x {0}. Therefore 

\z(A)\ $ e"**W £ 1 / d<?i...d<7 fc ££ e-^l^ e- A '^l Jj( J] ||T fl | 

fc ^ "A links j=l BeGj 

(6.8) 

The shortcut J dQ means a sum over the number m of transitions, a sum over transitions 
B\, . . . , B m , an integral over ordered times n, . . . , r m , and a sum over (m — 1) vertical 
segments that link {B{ x {r^}} together. 
We define 

z(G) = e~ A °M JJ ||Tb|| e^^+^l , (6.9) 
see 

where |£7| is the total length of the vertical segments of Q. If /3A/£ U 2^(log M + r), we 
can write 



E^pf/^/ d<W))*> (6.10) 



k>0 



where the first sum corresponds to the number of gathers linked to Q x {0}, and the 
second sum is the number of independent gathers. The shortcut j & dQ is identical to 
J dQ, except for the absence of an integral over n, which is set to 0; integrals over <S' and 
(3" are similar. 

One easily obtains an upper bound for the gathers with a unique transition: 

f dg~z{G)= V \\K B \\e 2utv ^ M+T ^ < (2-R + l)l.K"|| c (6.11) 

J& B,B30 

with c = 2vH. u (2R + l)"(log M + r); this is smaller than Ao if c is large enough in the 
assumptions of Theorem 4.4 . 

For general gathers, we proceed by induction. First, 

/ dGz(G) ^ (2R+l) v (\\K\\ c +\\K'\\ e +\\K"\\ c ). (6.12) 

Next, we use the recursion inequality: 

f dG~z(G) < £ || TB || e 2^(logM + r)|B| 
&m B,B30 

. / 

k yG_B ' 



E^( 2 E fdre- A - / (6.13) 
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Integrating over r, and since ||.K"|| c /Ao ^ 1 for a large enough c, we get 

dG~z{G) ^ \\T B \\e 2vr ^ M ^)\B\ e 2\B\ 

B,B90 

^ (2R + iy(\\K\\ c +\\K'\\ c +\\K"\\ c ). (6.14) 

This holds independently of m. This allows to estimate the integral over gathers that 
contain at least two transitions of type K or K" . Let & m C <3 m be gathers with at least 
one transition of type K or K' . One easily obtains 

/ 6Qz(Q) < (2R + 1) U (\\K\\ C +\\K"\\ C ). (6.15) 

Then the integral over gathers with two transitions of type K or K" can be done by 
integrating first on the time for such a transition, then over vertical segments and gathers 
at their ends, at least one of which must belong to ®^ n _ 1 . We obtain 

[ P dr V (\\K B \\ + \\KU)e 2 ^^ M+Tm f dr' e' A ° T ' [ dG~z{Q) 
Jo Jt 4 ' Jo J&L , 



B,BCQ 



ISO ,eB Jo J *"-™ 



< /?iei (l|Jf|k+ A l|g " lk>2 . (6.i6) 

Plugging these estimates in ( |6.10[) , one easily gets 

z(A)e T ^ ^ e-^o^e 3 ^. (6.17) 

Exponential decay of the weights of the contours is now clear. 

The bound on the derivative can be proven in the same way. Looking at ( |6.6D , we see 
that the integrand gets a factor bounded by (3\A\ sup x ^ UJ j \-^-^x{^u(x))\- 

6.3. Other properties of the weights. The weight of the contours can be viewed as a 
series in powers of {Kb}-, {K' b }, {K b }. Since it is absolutely convergent uniformly in K, 
K f , K" (provided they be small enough), we have by the dominated convergence theorem 

lim z(A) = 0. (6.18) 
k,k',k"-^o 

Analyticity of z{A) as function of fj, and (5 is clear, as well as a function of r/ if we add 
a new perturbation r/L, in a neighborhood of that depends on Periodicity is also 

obvious. 

Acknowledgments: It is a pleasure to thank Roberto Fernandez, Roman Kotecky, and 
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